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Abstract 

The clebsch potential approach to fluid lagrangians is developed in 
order to establish contact with other approaches to fluids. Three vari- 
ants of the perfect fluid approach are looked at. The first is an explicit 
linear lagrangian constructed directly from the clebsch potentials, this 
has fixed equation of state and explicit expression for the pressure but 
is less general than a perfect fluid. The second is lagrangians more 
general than that of a perfect fluid which are constructed from higher 
powers of the comoving vector. The third is lagrangians depending on 
two vector fields which can represent both density flow and entropy 
flow. 
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1 Introduction 



1.1 Motivation 

The motivation of this paper is to provided a unified approach to the various 
ways that fluids are described in physics. In particular the methods used 
by relativists, fluid mechanists, and nuclear physicists have grown distinct. 
In many areas of physics a unified approach is provided by the lagrange 
method, which for fluids is developed here. 

1.2 Methodology 

The methodolgy used is first to simplify a perfect fluid in order to investigate 
if methods of field theory can be applied to it; and then to generalize a perfect 
fluid to try and establish contact with more physical fluids. 

1.3 Other approaches to fluids 

A perfect fluid has a variational formulation [18, 36, 38, 39] which uses 
the first law of thermodynamics. In such a formulation clebsch potentials 
[10, 24, 4, 5, 17, 37, 13] for the comoving fluid vector field are used. Here 
this approach is both applied to less general fluids and to more general 
fluids. Other approaches to fluids include the following twelve. The first 
uses lagragians dependent on combinations of clebsh potentials which do not 
necessarily form a vector [34] . The second is that the comoving vector can 
be thought of as U a = x a , so that a perfect fluid is a type of generalization 
of a point particle, then there turns out to be a fluid generalization of a 
membrane [ ] . The third is that the charge subsitution d a — > d a + ieA a can 
be applied to fluids as well as fields and this leads to a model of symmetry 
breaking [32]. The fourth is that the navier-stokes equation has a lagrangian 
formulation [8, 11, 12, 15, 40], but the lagrangian has different measure and 
also image fields. The fifth is that hydrodynamics can be expressed using a 
grad expansion [16, 22, 21, 14] which needs an entropy vector. The sixth is 
that contemporary bjorken models use the grad expansion [7, 27, 28, 29, 20]. 
The seventh is fluid plasmas [3, 1]. The eighth is elastic models [2] §3, where 
the density rather than the pressure is used as the lagrangian. The ninth is 
other quantization methods such as brst and path integral applied to fluids 
j ]. The tenth is superfluids [ ]. The eleventh is spinning fluids [19, 23, 30]. 
The twelfth is cosmology [25], where clebsch potentials have been used. 
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1.4 Conventions 

The word potential is disambiguated by refereing to potentials for a vector 
field as clebsch potentials and potentials that occur in lagrange theory as 
coefficient functions. When a measure is suppressed it is f yf^^gdx 4 not J dr 
unless otherwise stated. \i is density and V is the pressure, q is a clebsch 
potential, a is used for a clebsch potential, a pauli matrix and the shear of 
a vector, to disambiguate the pauli matrix is always a p and the shear labels 

with which vector it is with respect to ex. Capital II indicates a momentum 
with respect to the proper time r not the coordinate time t. a,b,c, . . . are 
spacetime indices, i,j, k, . . . label sets of fields and momenta, and l,k,... 
label constraints. The signature is — h ++• 

2 The perfect fluid 

For a perfect fluid the lagrangian is taken to be the pressure C 
the action is 

I = J dx 4 V. 

The clebsch potentials are given by 

hV a = W a = a a + 6s a , V a V a = -1, (2) 

where if more potentials are needed it is straightforward to instate them; 
there are several sign conventions for 2. The clebsch potentials are sometimes 
given names: a is called the higgs because it has a similar role to the higgs 
field in symmetry breaking using fluids [32, 34], 6 is called the thermasy and 
s the entropy [38]. Variation is achieved via the first law of thermodynamics 

5V = n5h - nT5s = -nV a 5W a - nT5s, nh = ^ + V, (3) 

where n is the particle number and h is the enthalpy. Metrical variation 
yields the stress 

T ab = {fi + V)V a V b + Vg ab , (4) 
the nother currents j a = 5I/5q a are 

f a = -nV a , j% = 0, f s =-n9V a , (5) 

variation with respect to the clebsch potentials gives 

{nV a ) a = h + nG = 0, s = 0, 9 = T, (6) 



V, and 
(1) 



3 



where = V a a is the vectors expansion: thus the consevation of the nother 
currents 5 gives the same equations 6 as varying the clebsch potentials; the 
normalization condition V a V a = — 1 and 6 give 

a = -h. (7) 

The bianchi identity is 

T^ b = nW a + V' a , (8) 

subsituting for W using 2 and for V using 3 this vanishes identically. If 
one attempts to apply existing scalar field fourier oscillator quantization 
procedures to the above there is the equation 

W^ a = Ua+e a s a +eUs = {hV a ) a = h+hQ = h-^h = h(^n(^^ , (9) 

and if this vanishes the enthalpy n is proportional to the particle number n, 
for an example of the see the next section §3. The pressure V and density 
\x are only implicity defined in terms of the clebsch potentials so it is not 
clear what operators should correspond to them. Another possibility is to 
note that 6 are first order differential equations and to try and replace them 
with spinorial equations; however this would require a spinorial absolute 
derivative in place of the vectorial absolute derivative, see [31]§4.4. 
The canonical clebsch momenta are given by II 1 = 51 /5q l 

IF = -n, if = 0, IT = -710, (10) 

and these allow the nother currents 5 to be expressed as 

f q =n q v a . (ii) 

The standard poisson bracket is defined by 

r A n . 5 A 5B 5A5B 

(12) 

where i which labels each field is summed; the dirac matrix is defined by 

C« = {&,^}, (13) 
and the dirac bracket is defined by 

{A, B} * = {A, B} - {A, fa} Inv (C tK ) B} , (14) 
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where Inv(C tK ) denotes the inverse of C M . For a perfect fluid the constraints 
are 

0! = it - #rr, (h = if, (15) 

and the dirac matrix 13 becomes 

C LK = -ia p2 n°, InvC tK = +i^r, (16) 
where a P 2 is the pauli matrix 

^=(1 J)' " p2 "(+* o')' CTp3= (o (17) 

The dirac bracket 14 for 15 is 

acting on the clebsch momenta and clebsch potentials 

{IF , a}* = -1, {IT 9 , #}* = 0, {IT, s}* = -1, 

n CT {a, 0}* = -e, ir{0,s}* = -i, n CT {CT,s}* = o, (19) 

replacing the dirac brackets by quantum commutators there does not seem 
to be an easy representation of the resulting algebra as discussed in [33] . To 
quantize one replaces the dirac bracket with a quantum commutator 

{A, B} * -> — \AB - bA] , (20) 
in I J 

in the present case the second equation of 19 suggests that U s = 0, which is 
assumed from now on. There is the problem of how to realize the coordinate 
commutation relations in the last three of 19 which does not seem possible 
using the pauli matrices 17 and so this is left for now. Suppose that IF is 
replaced by a differential operator 

IT = -ihd x , (21) 

possibilities for x include spacetime coordinates x a , the particle number n 
or the proper time r which restrict a as to a lesser extent does the dirac 
operator ^y a d a so the simplest choice is taken d a . To proceed it is necessary 
to have an explicit hamiltonian which does not exist in the present case as 
the pressure V and the density \i are not directly expressible in terms of the 
clebsch potentials. 
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3 The explicit linear lagrangian 

The lagrangian is taken to be linear in W 2 

£ = -^(n,q)W 2 -U(n,q), (22) 

where T{n, q) and U(n, q) are the first and zeroth order coefficient functions 
of the particle number and clebsch potentials respectively. Metric variation 
yields 

T ab = TW a W b + Cg ab = Th 2 V a V b + Cg ab , (23) 
requiring this stress to be that of a perfect fluid places the restriction 

fh 2 = nh = T + n, or T = ^, (24) 

again requiring that the pressure is the lagrangian and using W 2 = — h 2 
gives 

C = V=\nh-l3= l -{V + l j)-l3, (25) 

which yields the linear equation of state 

V = fi-2U, (26) 

this is a very restrictive equation of state as it does not have even simple 
cases such as V = (7 — as examples. Variation with respect to the 
clebsch potentials gives 

n + n9-O iCT = 0, ns + O >e = 0, n9 - 0, s + 0O I<7 = 0, (27) 

where 

% = V, q ~ \h 2 F q , (28) 

and 24 has been used. Using the equations of motion 27 the bianchi identities 
8 become 

T% = \h 2 T a - (w + \h 2 ^)j s a - n9 - a a l3 a - 8 a U e , (29) 

for the thermodynamical case 6 and 27 give zeroth order coefficient function 
15 = nTs and then the bianchi identity reduces to T a = or n = kh, 
thus what in lagrangian 22 looked like an arbitrary function T has been 
reduced to a constant and this will be further considered in the next section. 
Variation with respect to everything else except the particle number n gives 
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the same as the implicit approach of the previous section §2. For variation 
with respect to the particle number n there are three choices: the first is 
to simply do it in which case, for lagrangians in which n is separated and 
linear, the vectors W, V are forced to be null and the system is no longer 
that of a perfect fluid, the second is to ignore variation with respect to n, for 
lagrangians in which n is separated and linear this amounts to assuming a 
first law of the form 3 and one is essentially working in the implicit formalism 
of the previous section §2, the third is to alter lagrangians with coefficient 
functions F, 15 that depend on n, then variation with respect to n can be 
chosen so that h 2 T. n = — 2l5 :Tl . Using 10 for n, 7 for h and 2 for V 2 the 
lagrangian is 

£ = V = ^all a - U, (30) 

using 27 the hamiltonian is 

H= -&U a + U- -U a U e , (31) 
2 n 

when I5 t $ = the hamiltonian equals the density Ti = fi. For this explicit 
form 31 it is possible to use the operator substitution 21 so that 

W& = -^iha^ a + m = 0, (32) 

taking a^f a = 32 becomes 

ihi> = 2m, (33) 

integrating 

* = Aexp J Udr^j , (34) 

this wavefunction turns out to be too restrictive to be of any use. 

4 Simplest three clebsch potential fluid 

The explicit linear thermodynamic lagrangian is 

C = -\wl-nTs, (35) 
see the remarks after 29. Varying with respect to the potentials gives 
Wt = Ua + 9 a s a + ens = 0, s a (a a + 6s a ) = 0, 6 a (o- a + 6s a ) = hT. (36) 
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If these equations are thought of as functions of one variable then the second 
equation gives either s a = or W a = both of which are of no practical use: 
therefore two variables are used specifically to seek plane wave solutions. A 
simple choice is 

a = Aexpi (hot + k\x) , = i?exp at (hot + k±x) , s = C exp(l— a)i (k^t + k\x) , 

(37) 

which gives 

W a = i[k o ,k 1 ,0,0])C, K=(j^ + l-c^9s, (38) 

the first and second equations of motion 36 give 

{kl - k\)K = 0, (39) 

and the third gives 

{kl-kl)K=T- el (40) 

thus the system is decoupled, in other words as 39 forces it to be null the 
terms in the first equation of 36 vanish separately, and by 40 the system is 
at zero temperature. Working through in the same manner with the clebsch 
potentials having both left and right movers, i.e. like 48, the same problems 
arise and there is no indication that mixed movers can generate non-zero 
temperature. Working through without the restriction that the clebsch po- 
tentials are co-directional produces equations too general to proceed with. 
The hamiltonian is 

H = -^&W -TsU a , (41) 
in the present case a = h = n = —W so that 

n = ^rr - Tsir (42) 

using the substitution 21 with x = a gives 

ih 

^ = __* ctct + Ts ^ = 0, (43) 

which has solution 

* = Aa exp i- 2 ^) + B, (44) 



h J 



which is again too restrictive. 
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5 One clebsch potential fluid with two coefficient 
functions 



The lagrangian is taken to be 

C = -^(a)a 2 a -U(a), (45) 

the metric stress is 

T ab = To a o h + £g a b, (46) 

V, /x, n, h are recovered in the same way as in the last section. Varying with 
respect to a 

TUo + T a a a + I3 a = 0, (47) 

which for simple T has simple spherically symmetric solutions; for wave 
solutions use zeroth order coefficient function 13 = m 2 a 2 /2 and 

a± = A + exp(ik ■ x) ± A_ exp(—ik ■ x), (48) 

with a = 0"+, 47 becomes 

[Tkl + m 2 )a + T'k 2 a o 2 _ = 0, (49) 

the last term forces either or A- to vanish. 

6 One fluid described using one vector field 

In the one fluid one vector approach one hopes to find a lagrangian which 
recovers as much of the stress [2G]eq.22.16d as possible 

Tab = {^ + V -Ce)U a U b -2 V aab+2q {a U b) + {V -^e)g ab , 
S a = nsUa + j;, (50) 

where \x is the density, V is the pressure, £ > is the coefficient of bulk 

viscosity, rj > is the coefficient of dynamic viscosity, a is the shear, q 
is the heat flux, S a and s are the entropy vector and scalar and T is the 
temperature. Consider a fluid lagrangian dependent on one velocity which 
can be expanded 

C = C(V) = C(V 0+ ,V 1 ,V 1+ ,V 2 ,...) (51) 
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where 

C(V 0+ ) = k 0+ V, (52) 

ay 1 ) = k l Q = k l v% 

c(v 1+ ) = k 1+ re, 

cy 2 ) = kfe + klR ah v a v b + kluj 2 + kja 2 + fcf e 2 + k\v* a 



so that the integer superscript on the velocity V indicates the power in 
which it occurs in the lagrangian, the meaning of the superscript + will 
become apparent later. Here just the first three terms are considered, terms 
of second or Raychadhuri order and higher are ignored, as are any auxiliary, 
image, entropy or electromagnetic fields. From a technical point of view just 
the first term in this expansion has been considered in the previous section 
§2 and in this section the next two terms are considered. Metrical variation 
gives the stress 

T ab = (n + V) {k 0+ + k 1+ @) V a V b - 2*;%.,,) + Cg ab , (53) 

53 does not bare much resemblance to 50 except that there is a term similar 
to bulk viscosity appearing. Variation with respect to the clebsch potentials, 
particularly of Cy 1 ), gives long expressions which are not helpful. Variation 
with respect to the clebsch velocities gives the momenta 

II* = -n (k 0+ + k 1+ @) , IT = m*, if = 0, (54) 

and these give two constraints 

<p 1= U s -eU^, 2 = n e . (55) 

In the present case, as the constraints 55 are the same as for the perfect 
fluid, the dirac bracket between the coordinates and momenta and thus the 
quantum relations are the same as for the perfect fluid as given above and 
in [33]. 

7 One fluid described using two vector fields 

If one tries to incorporate entropy by choosing it to move in the same di- 
rection as fluid flow then the entropy is proportional to the enthalpy and 
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s = kh and the first law becomes 

Sp = nbh - nTbs = -nV a 5W a - nTkbh 

= -nV a 5W a + nTkV a 5W a = -n*V a bW a , 

where n* = (1 — kT)n, (56) 

so the problem becomes the same as before with the particle number n 
replaced by n* . In the one fluid two vectors approach one proceeds as before 
except the clebcsh decomposition 2 is replaced by 

hV a = W a = a a + nxa, V a V a = -1, 

lU a = X a = a a + p-ya, U a U a = -l, (57) 

and the first law 3 is replaced by 

5 P = nbh - nT5S = -nV a 5W a + nTU a 5X a . (58) 

Metrical variation gives 

T ab = nhV a V b - nT£U a U b + Vg ab , (59) 

when either T or £ —> the second term vanishes. 59 can be re-written as 

T ab = {V + n) (v a V b - Yi XaX ^) + V ^b, (60) 

Two choices firstly the free case where the clebsch potential are all inde- 
pendent, secondly the thermodynamic case where the clebsh potential are 
chosen so that the thermasy and entropy change are non- vanishing. The 
free case essentially duplicates the one vector case of the perfect fluid §2. 
For the thermodynamic case choose 

e = v = (3, s = x = l, (61) 
in 57. Defining two absolute derivatives 

T a bc... = V e T a bc...;e, T' abc = U e T a b c ...-e, (62) 

variation with respect to a, a, 6 and s gives 

n + n 6= 0, (nT)' + nT Q= 0, s = Ts', 6 = TO' , (63) 

respectively. The relationship between dot and dash derivatives is needed, 
note 

X a = (a- a) a + W a , (64) 
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thus for an arbitrary tensor Aiftc... there is the relationship between the 
derivatives 

£A' abc _ = (a- a) e A abc ... e + hA abc ..., (65) 
using 65, 63 becomes 

hQ-£Q= (ln(nT)) a (a - a) a + h (ln(T))° , 

. = t(a - a) a sg ■ = £(a - a) a e a 

S i-Th ' " l-Th ' 1 ' 

In the limit that the two vectors coincide a — ► a and 64 and 66 give T = 
s = 6 = so that the thermasy relation is not recovered unless T = 0. The 
dot momenta are the same as that given by 10, the dashed momenta are 

IP (dashed) = nT, II s (dashed) = n6T, (67) 

it is necessary to convert these to dot momenta using 

nl " w = = w ul{dashed) " ^^ n ^ dashed )' ^ 68 ) 

where i is not summed and / is an undetermined function of the clebsch 
potentials; collecting together 10,67,68 gives the total momenta 

IT = -n, IT* = nTf, U e = 0, IT = -n6 + n9Tf, (69) 

and the three constraints 

X = n s - ew - 9n a , <f> 2 = n 9 , 3 = n a + Tfu a , (70) 

the dirac matrix is 

c 12 = -w-u a ,c 13 = (- 6 ^ + e^ + e^jTW,c 2 3 = -^Tu^ (71) 

the inverse InvC can be found, but there is no simplification in further 
expressions. 



8 Conclusion 

In §2 the lagrangian approach to perfect fluids was presented, the clebsch 
potentials contain more information than is needed to describe the system 
so that it is constrained, dirac constraint analysis removes the superfluous 
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degrees of freedom resulting in the algebra 19, see also [33], which when 
quantized via 20 does not seem to lead to a recognizable algebra; also as the 
pressure and density are implicit rather than explicit functions of the clebsch 
potentials it is not clear what quantum operators should represent them. In 
§3 to overcome lack of explicit expressions for the pressure and density an 
explicit lagrangian 22 was studied, this has a very restrictive equation of 
state 26, and for thermodynamic choices of the functions reduces to the 
example of the next section. In §4 the simplest three potential lagrangian 
35 was studied, it seems to lead to an unrealistic quantum theory 43. In §5 
a one potential lagrangian 43 was studied, it is just a simple generalization 
of the klein gordon lagrangian, but the generalization is obstructive enough 
to prevent successful investigation of it by the fourier oscillator method. In 
§6 the perfect fluid is generalized so that it depends on higher powers of 
the comoving fluid, a term similar to the bulk viscosity appears and the 
momentum constraints are the same as those for a perfect fluid. In §7 the 
perfect fluid is generalized to include two vector fields in the hope that these 
can represent both density and entropy flow, equating some of the clebsch 
potentials in each vector leads to plausible thermodynamic equations 66. 

9 Acknowledgements 

I would like to thank Tom Bridges, Alex Craik and Mattias Marklund for 
their interest in parts of this work. 

References 

[1] A. Achterberg, Variational principle for relativistic magnetohydrody- 
namics, P%s.-Rev.428(1983)2449-2458. 2 

[2] Lars Andersson, Robert Bieg and Bernd G. Schmidt, Static 
Self-Gravitating Elastic Bodies in Einstein Gravity (2008), 
Math.Rev. 2009f:83034, Comm.Pure Appl. Math. 61(2008)988-1023. 2 

[3] Angelo Marcello Anile, Relativistic Fluids and magneto-fluids, Cam- 
bridge University Press (1989). 2 

[4] P.R. Baldwin, (1993) Constructing Clebsch Potentials for Vector Fields. 
Print Physics/ Akron. 2 

[5] Eugene Balkovsky, Some notes on the Clebsch representation for incom- 
pressible fluids. Math.Rev. 94k: 76026 Phys. Lett. A 186(1994)135-136. 2 



13 



Arjun Berera & David Hochberg, Gauge fixing, BRS invariance and 
Ward identities for randomly stirred flows, Nucl. Phys. 5 814(2009)522- 
548. 2 

James D. Bjorken, Highly relativistic nucleus-nucleus collisions: The 
central rapidity region, Phys. Rev. D 27(1983)140. 2 

Carlos Cartes, Miguel D. Bustamante & Marc Etienne Brachet, Gener- 
alized Euler-Lagrangian description of Navier-Stokes dynamics, Physics 
of Fluids 19(2007)007101. 2 

Brandon Carter & I.M. Khalatnikov Momentum, vorticity, and helic- 
ity in covariant superfluid dynamics. Math. Rev. 93m:83046 Ann. Physics 
219(1992)243-265. 2 

Alfred Clebsch, Uber eine allgemeine Transformation der Hydrody- 
namischen Gleichungen, Crelle's Journal, liv (1859). Reine Angnew. 
Math. ("Crelle")56(1859)l. 2 

Peter Constantin, An Eulerian-Lagrangian Approach 
to the Navier-Stokes Equations, Math.Rev.2002m:76023 
Comm. Math. Phys.216(2001), 663-686. 2 

Hugh L. Dryden, Francis P. Murnaghan and H. Bateman, (1956) Hy- 
drodynamics, Dover, Math.Rev. 0077307 (17,1019f). 2 

Carl Eckart, Variational Principles of Hydrodynamics, 
Math.Rev.0114428 (22#5249), Physics of Fluids 3(1960)421 
esp. Appendix. 2 

A. El, Z. Xu & C. Greiner, Third-order relativistic dissipative hydro- 
dynamics, 0907.4500 2 

Bruce A. Finlayson, Existence of variational principles for the 
Navier-Stokes equation, Math.Rev.0326189 (48#4534), Phys. Fluids 
15(1977)963-967. 2 

Harold Grad, On the kinetic theory of rarefied gases, Math.Rev. 0033674 
(11,4739), Commun.Pure Appl.Math.2(194:9)331. 2 

C. Robin Graham and Frank S. Henyey, Clebsch representation near 
points where the vorticity vanishes, Math. Rev. 2000k: 76033 Physics of 
Fluids 12(2000)744-746 2 



14 



[18] R. Hargreaves, A Pressure-integral as Kinetic Potential, 
PM.Met 5 .16(1908)436-444. 2 

[19] Friedrich W. Hehl, Paul von der Heyde, G. David Kerlich and James M. 
Nester, Math.Rev. 0439001 (55#11902), Rev.Mod.Phys.48(1976)393. 2 

[20] Michal P. Heller, Romuald A. Janik and R. Peschanski, Hydrodynamic 
Flow of the Quark-Gluon Plasma and Gauge/Gravity Correspondence, 
0811.3113 v2 2 

[21] William A. Hiscock and Lee Lindblom, Math.Rev.88e:83056, 
Phy -s. Rev. .035(1987)3723-3737. 2 

[22] W. Israel and J.M. Stewart, Transient relativistic thermodynamics and 
kinetic theory, Math.Rev. 80d:8001 1, Ann.Phys.N. 7.118(1979)341. 2 

[23] Roman Jackiw, Description of Vorticity by Grassmann Variables 
and an Extension to Supersymmetry, Fluctuating paths and fields, 
553-563, World Sci.Pub., River Edge NJ (2001), Math.Rev. 1868089 
physics/0010079 v2 2 

[24] Horace Lamb, Hydrodynamics,p.248, CUP, Cambridge 1932, 
Math.Rev.96£-7600L 2 

[25] V.G. Lapchiniskii and V.A. Rubakov, T/teor.Mai/t.P%s.33(1977)1076. 
2 

[26] Charles W. Misner, Kip S. Thorne & John Archibald Wheeler, Gravita- 
tion. W.H.Freeman & Company, New York (1970), Math.Rev. 0418833 
(54#6869). 9 

[27] Ian G. Moss, Superfluidity in Super- Yang-Mills Theory, 0712.2568 2 

[28] Azwinndini Muronga, Causal theories of dissipative relativistic fluid 
dynamics for nuclear collisions, Phys.Rev. C 69(2004)034903. 2 

[29] Shin Nakamura and Sang-Jin Sin, A holographic dual of hydrodynam- 
ics, Math.Rev.2007e:81i20, JEEP 09(2006)020. 2 

[30] Tino S. Nyawelo, Jan-Willem van Holten and Stefan Groot Nibbe- 
link, Relativistic fluid mechanics, Kahler manifolds and supersymmetry, 
hep-th/0307283, Math.Rev.2005b:83053 2 

[31] Roger Penrose and Wolfgang Rindler, Spinor calculus and relativistic 
fields, CUP (1987), ISBN 0521-33707-0 4 



15 



[32] Mark D.Roberts, Fluid Symmetry Breaking II: Velocity Potential 
Method, Hadronic J.20(1997)73-84, hep-th/9904079 2, 3 

[33] Mark D. Roberts, The Quantum Commutator Algebra of a Perfect 
Fluid. Math. Rev. 2000j:81120 Mathematical Physics, Analysis and Ge- 
ometry 1(1999)367-373. gr-qc/9810089 5, 10, 13 

[34] Mark D. Roberts, A Generalized Higgs Model. Physics Es- 
says(2006) September Vol.19, No. 3, p.1-5. hep-th/9904080 2, 3 

[35] Mark D. Roberts, A Fluid Generalization of Membranes, 
hep-th/0406164 2 

[36] Hanno Rund, D.R. Wells and L.C. Hawkins, J. Plasma 
Phys. 20(1978)329-344. 2 

[37] Hanno Rund, Clebsch Representations and Relativistic Dynamical Sys- 
tems, Math.Rev.80d:83054, Arch. Rat. Mech. Anal. 71 (1979) 199-220. 2 

[38] R.L. Salinger and Gerald B. Whitham, Variational Principles in Con- 
tinuum Mechanics, Proc. Roy.Soc. London A305(1960)l-25. 2, 3 

[39] Abraham Haskel Taub, General Relativistic Variational Principle for 
Perfect Fluids, Math.Rev. 0063824 (16, 185b) Phys.Rev.6(1954)U68- 
1470. 2 

[40] J.R. Usher & Alex D.D. Craik, Nonlinear wave interactions in shear 
flows. Part 1. A Variational formulation, J. Fluid Mech. 66 (1974) 209- 
221. 2 

[41] Victor Yakhot and Vladimir Zakharov, Hidden conservation laws in 
hydrodynamics, Math.Rev. 94a:76041, Physica D 64(1993)379-394. 



16 



